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Abstract. We consider the effective dynamics obtained by double-passing a far- 
detuned laser probe through a large atomic spin system. The net result of the atom- 
field interaction is a type of coherent positive feedback that amplifies the values of 
selected spin observables. An effective equation of motion for the atomic system is 
presented, and an approximate 2-parameter model of the dynamics is developed that 
should provide a viable approach to modeling even the extremely large spin systems, 
with N ^> 1 atoms, encountered under typical laboratory conditions. 

Combining the nonlinear dynamics that result from the positive feedback with 
continuous observation of the atomic spin offers an improvement in quantum parameter 
estimation. We explore the possibility of reaching the Heisenberg uncertainty scaling 
in atomic magnetometry without the need for any appreciable spin-squeezing by 
analyzing our system via the quantum Cramer-Rao inequality. Finally, we develop 
a realistic quantum parameter estimator for atomic magnetometry that is based on 
a two-parameter family of Gaussian states and investigate the performance of this 
estimator through numerical simulations. In doing so, we identify several issues, such 
as numerical convergence and the reudction of estimator bias, that must be addressed 
when incorporating our parameter estimation methods into an actual laboratory 
setting. 
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1. Introduction 

Real-time feedback plays an essential role in controlling the behavior of physical systems 
that must respond to unanticipated events on the same time-scale as their own dynamics 
[H [2], |3j . This basic engineering premise shows up everywhere: mechanical actuation 
based on situational detection provides the basis for devices like the air handling system 
in a typical optics lab as well as the autopilot on a typical airliner; electrical feedback lies 
at the heart of error correction in communication systems; and biochemical feedback 
implements metabolic control and maintains homeostasis in just about every living 
organism. Another example, and one which plays an important role in this work, involves 
the amplification of signals using positive feedback [J]. In the case of amplification, 
feedback must be implemented in real-time because a bona fide amplifying device must 
operate on signals whose form and value is not known in advance. 

Engineering quantum mechanical systems shares much in common with its classical 
counterpart in that optimal [5j E] and real-time feedback also provides an essential tool 
for controlling behavior in the face of unpredictable fluctuations. Examples include 
stabilizing Heiseneberg uncertainty to achieve deterministic state preparation [7J [8], 
implementing optimal measurements for quantum state discrimination [9], quantum 
error correction [10] , and affecting the motion [TH [12] or internal states [13] of single 
particles or qubits, to name a few. Quantum feedback does differ from its classical 
analog, however. In the classical picture of feedback, one continually measures the 
system to determine its deviation from the intended behavior by constructing an error 
signal and then manipulates the system to put it back on track by minimizing the 
error. Every part of that process is governed by the same type of physical laws, so 
to speak; the measurements, actuation, external influences, etc., can all be modeled 
using classical mechanics, thermodynamics, information theory, etc. The situation is 
different in quantum mechanics, where the system undergoes substantially different 
types of dynamics depending on whether it is being observed or not, which is significant 
because it implies that there are qualitatively different ways to enact feedback. Quantum 
feedback control is therefore typically broken down into different forms: 

(i) Measurement- based feedback: continuous observation of the quantum system is 
performed by entangling it with an auxiliary meter, such as the electromagnetic 
field, and measuring the meter continuously in time to gain information about the 
target system [HI [7]. A control Hamiltonian is then modulated by changing the 
values of classical parameters (such as field strengths) in the Hamiltonian as a time- 
dependent function based on the evolving measurement record [15J. The feedback 
controller is implemented using classical signal processing techniques to interpret 
the continuous measurement record and determine the appropriate feedback signal 
from a classical control law. Feedback is applied by modulating the system's 
quantum Hamiltonian function of time. 

(ii) Coherent feedback: the quantum system of interest is coupled coherently to an 
auxiliary quantum system according to some Hamiltonian. The auxiliary system 
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directly plays the role of the feedback controller [THl HZ]. Its dynamical state 
changes based upon that of the target system as a result of the coherent interaction 
that couples them, and the state of the system likewise depends upon that of the 
controller (either through the same interaction, or some additional coupling). In 
this sense, the dynamics of the target system are modified based on its own state 
via the external coupling. The coupling Hamiltonian responsible for the feedback 
is generally not modified clS db function of time. 

In this paper, we consider a hybrid type of control process that combines aspects from 
both of the two varieties of feedback just described. A laser field that is double-passed 
through an atomic spin system (q.v., Figure CD) [HI [19j [20], ED, E2] provides a type 
of coherent feedback interaction. On the first pass of the light through the atoms, 
the polarization state of the field is modified in a manner that depends upon the 
atomic spin. Then on the second pass, the interaction is configured such that the 
light drives the atoms in a way that depends upon the optical polarization. The net 
result is a positive-feedback amplification of selected atomic spin observables. Both 
of the atom-field interactions are coherent in nature, but the field is also an infinite 
reservoir that would decohere the atomic sample if otherwise neglected. Combining 
the atom-field interaction with continuous measurement of the twice-scattered field, 
however, helps to limit this decoherence. Measurement also enables the estimation of 
atomic spin observables, which in turn can be used to estimate the strength of external 
influences on the atomic spin. To demonstrate the utility of this system for quantum 
parameter estimation, we consider the problem of determining the magnitude of an 
external magnetic field. The field estimation procedure is designed in such a way that 
the positive coherent feedback caused by the optical field amplifies the spin dynamics 
that result from the magnetic field in a manner that improves the overall sensitivity 
[2H122]. 

The remainder of the paper is organized as follows. Section (T5]) develops an equation 
of motion for the coherent atom-field feedback interaction in the form of a stochastic 
time-evolution operator for the joint system. Section ([3]) develops the so-called quantum 
filter that allows one to estimate atomic observables based on continuous detection of 
the scattered field. Section (@J utilizes the filtering equations to develop an estimator for 
determining the strength of an applied magnetic field based on the induced dynamics 
of the atomic system in the combined presence of the magnetic and coherent feedback 
interactions. 

2. The Physical System 

Consider a physical system comprised of a cloud of iV atoms coupled to the 
electromagnetic field, as depicted by the schematic in Fig. ([1]). The degrees of freedom of 
interest are provided by the hyperfine angular momentum of the atomic system, which 
arises from the coupling of the spin and orbital angular momentum of the atom's valence 
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Figure 1. (color online) The collective spin angular momentum of N atomic spins is 
coupled to a far-detuned laser field that is double-passed through the atomic sample. 
On the first pass through the atoms, the field is linearly polarized and therefore acquires 
a Faraday rotation proportional to the z-component of the atomic spin. On the second 
pass, the Faraday rotation is transformed into optical helicity, which the atoms perceive 
as a fictitious magnetic field. Symbol defintitions: (PBS) polarizing beam splitter, 
(PD) photodetector, (A/2) half-waveplate, (A/4) quarter- waveplate. 



electrons with its nuclear spin 

f(j) = g(i) + i(j) + i(j). (1) 

Here, the superscript j = 1, . . . , N indicates on which of the N atoms the corresponding 
operator acts. A basis for the atomic Hilbert space J%a is then provided by the 
simultaneous eigenkets of the total spin and its z-projection 

(.f 'f ./}./„,) IrfSfj I 1 )./}•///,) (2) 

fz ] \fji m j) = ^ m j\fj, m j) (3) 

for each of the j = 1, . . . , N atoms. In many laboratory situations, however, the 
dynamic quantum state of the atoms resides in a small sub-space of J#a called the 
completely symmetric representation. Such states are obtained when the total spin 
quantum number is the same for all of the atoms fj = fj> = f and the state of the N 
atoms remains invariant under a permutation of particle labels Vjj>[p] = p [231 121] • The 
effective atomic Hilbert space is then spanned by the eigenkets 

F Z \F,M) = hM\F,M), M=-F,...,+F (4) 

of the collective spin operators 

N 

F, = J2ff ] - (5) 

Dynamics that are generated by functions of the collective spin operators preserve 
the total spin quantum number, and for spin-polarized atomic systems such as those 
considered here, its value corresponds to the maximal angular momentum states 
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F = Nf. The dimension of the symmetric atomic Hilbert space is therefore linear 
in the number of atoms gp — 2Nf + 1, and the atomic state will remain confined to this 
sub-space provided that the dynamics do not distinguish between the different atoms 
[23]. In this paper, only states and dynamical models that preserve symmetric collective 
states to good approximation are considered [23]. In practice, realizing such a model in a 
laboratory setting would require that the electromagnetic fields used to manipulate the 
atomic sample couple identically to each atom and that processes such as spontaneous 
emission, which do no preserve the completely symmetric representation, are heavily 
suppressed [23J. 

2.1. Atom-Field Coupling 

As also depicted in Fig. (JTJ, the atomic system is illuminated by a relatively intense, 
far-detuned laser. The traveling-wave electromagnetic field propagates initially along 
the atomic z-axis and is x-polarized, i.e., it is linearly polarized such that its electric 
field oscillates along the atomic x-axis. After its first interaction with the atomic system, 
the beam is steered such that it passes through the atomic system a second time, along 
the atomic y-axis [HI [19], [20]. Prior to this second interaction, the field is passed 
through polarization optics that are configured in the following manner: the x-polarized 
component of the field remains unaffected, while the ^-polarization component is phase- 
retarded by a quarter wavelength. Such a transformation is achieved in the laboratory 
by placing a A/4-plate in the beam path with its principal axis oriented along the x-axis. 
As a result of this configuration, any rotation of the laser polarization during its first 
pass through the atoms is converted into optical helicity on the second pass. The atoms 
perceive this optical helicity as a fictitious magnetic field, which drives rotations of the 
collective atomic angular momentum. 

Modeling the dynamics of the double-pass atom-field system requires that we 
develop an equation of motion for the joint state of the atoms and the field, for example 
by computing the time-evolution operator U t that brings an initial state po to the evolved 
state p t = UtPoU t at time t. To do so, we proceed according to the following steps: 

(i) a separate interaction Hamiltonian for each of the two passes of the far-detuned 
laser field through the atomic system is developed; 

(ii) a weak coupling limit is taken to obtain a Markov equation of motion generated by 
each of the separate interaction Hamiltonians; 

(iii) the two Markov-limit propagators are combined into a single coarse-grained Markov 
limit that describes the dynamics on time-scales that are slow compared to the time 
required for the laser field to propagate twice through the atomic system. 

Our intention here is not to present a detailed derivation of the stochastic time-evolution 
operator and quantum filtering equations that result from the procedure just described, 
but rather to demonstrate the capabilities of the measurement configuration in Fig. ([I]). 
For details on the derivations, we refer the interested reader to Ref. [22J. 
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2.2. The Hamiltonian for the Atom-Field System 

When the collective spin angular momentum of a multilevel atomic system interacts 
dispersively with a traveling wave laser field with wavevector k, the atomic spin couples 
to the two polarization modes of the electromagnetic field transverse to k according to 
the atomic-polarizability Hamiltonian 



N 



HA - iT/2 



Efo.t) (6) 



where rj is the position of atom j, is the dipole lowering operator for atom j, and 
A = cjl — k>A is the detuning of the laser from the relevant atomic resonance wa and F is 
the natural linewidth of that resonance. The atomic polarizability Hamiltonian can be 
obtained via second-order perturbation theory in which the electronic excited states of 
the atoms are adiabatically eliminated to obtain an effective Hamiltonian on its ground 
state ES]. 

To work with Eq. (jfJJ), it is convenient to treat two polarization modes orthogonal to 
the direction of laser propagation as a Schwinger-Bose field. When quantized in terms 
of a plane-wave mode decomposition, one obtains the Stokes operators: 



1 



s 



I 

= +2 ( fi y,<A,u, - a x lt A, J ( 8 ) 

1 (^L^x^ + a x>y J ( 9 ) 



s = — 



+ 2 ( a y,<A,^ - a x > y ,J , (10) 

expressed in terms of the Schrodinger-picture field annihilation operators, a X£J and a yu! , 
for the plane-wave modes with frequency uj and linear polarization along the x- and y- 
axes, respectively. But, as is evident from Eq. ([6]), the interaction Hamiltonian couples 
the atoms to the electric field as a function of time, suggesting that it would be more 
convenient to transform the Stokes operators from a plane-wave mode decomposition of 
the electromagetic field to operators that are labeled by time. Toward this end, and after 
moving into an appropriate interaction picture, it is convenient to define the following 
time-domain Stokes lowering operator 

st = -J g^a^a^-^dwi (11) 

where g(u) is an atom-field coupling parameter that is a function of the field amplitude, 
detuning, and the atomic transition dipole matrix element [7]. The Stokes operators 
corresponding to the first-pass and second-pass directions can then be expressed in terms 
of §t as follows: 

s Z;t = i (s\ -s t J and s y>t = - + s{) » ( 12 ) 
making it possible to view the polarization operators as field quadratures. 
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It can be shown that the interaction Hamiltonians for each pass of the probe light 
through the sample take the form [27] 

= +ihX 1 F z (13) 

HP = - h\ 2 F y (s t + 4) , (14) 

respectively. The coupling constants Ai and A2 have been introduced to allow one to 
scale the strength of interaction, and two distinct constants are introduced to account for 
the possibility that the first and second-pass atom field interaction strengths differ. In 
the simplest experimental setting, the laser field will have very nearly the same intensity 
and detuning on both passes, such that Ai = A2. 

Note that in getting from Eq. (jBj) to Eqs. (11311141) . it was assumed that the traveling- 
wave nature of the field allows one to average over the atomic position and that rank-two 
spherical tensor contributions to the dyad dcr can be neglected [281 [29]. I n practice, such 
an approximation depends heavily upon the choice of atomic species, laser detunings, 
and intensities [28, 29J. As spatial field modes other than those corresponding to the 
direction of the laser propagation have been omitted from the atom-field Hamiltonian, 
atomic spontaneous emission into non-paraxial modes has implicitly been neglected. 
Such an approximation can be justified if the laser is sufficiently far detuned, although 
in practice, care must be taken to ensure that such conditions are met J]] 



2.3. The Stochastic Propagator 

The objective now is to determine the time-evolution of the atom-field system subject 
to the two Hamiltonians, Eqs. (11311141) . which can be addressed by computing the time 
evolution operator on the joint atom- field system U t . Toward this end, an equation of 
motion for the propagator is developed by analyzing the two interaction Hamiltonians 
separately, and then combining them. Specifically, this involves calculating the 
increments dU^ and dU^ such that the time-evolution operator just after the first 
pass of the light through the atomic sample is given by 

U t+dt = U t + dUi 1) (15) 

and the time-evolution operator following the second pass through the atoms is 

U t+2dt = U t+dt + dlli 2) . (16) 

Then, an effective increment dU t is obtained that combines the actions of dUf and 

" (2) 

dUf' into a single update equation for Ut- 

Beginning by analyzing the first pass interaction, the propagator satisfies the 
equation of motion 

^ = -^m t = XiF z (s t -r s \)u t , (17) 

I Neglecting spontaneous emission is justified by the fact that we only consider the case of extremely 
weak measurements in which the atom-field coupling strength is not large enough to induce appreciable 
backaction on the atomic system by measuring the forward scattered laser field. 
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which can be formally integrated to give 

Ut = 1 - Ai / drP z (s T - 4) U T , (18) 
Jo 

from which the increment 5U^ = U t +s — U t can be defined 

ft+5 rt 

SU^ 1} = X 1 / drP z (s T - 4) U T -X 1 / drP z (s T - 4) U T (19) 
and used to obtain a differential equation for U t by taking the limit 6 — > 0. 

2.3.1. The van Hove Limit. The 5 — > limit in Eq. (JT9l) is evaluated by taking the 
atoms to be weakly coupled to the field, which one would expect to be the case for 
a finite coupling strength applied over an extremely short period of time St. This so- 
called "weak-coupling limit" is achieved by scaling the perturbation parameter A — > 0, 
however, several technical complications arise in doing so that require one to adopt 
the methods of quantum stochastic calculus [30l EU [7j. It is, of course, insufficient to 
take A —>■ simply on its own, as doing so would just turn off the interaction. The 
objective of the weak-coupling limit should be viewed in the sense of making a Markov 
approximation — it is to determine the effective dynamics of the system on time-scales 
that are long compared to the correlation time of the field fluctuations. 

In the theory of quantum stochastic calculus, such a limit is obtained by taking 
t — ► oo at the same time as A — > by rescaling time according to t i— > t/X 2 . This 
procedure, often called the van Hove limit, has the desired effect of exposing the 
effective dynamics of the system that result from its average interaction with the rapid 
fluctuations of the field. The fact that time should be rescaled as t/X 2 , and not some 
other function of A, is a byproduct of the correlation function for operators on the field 

[7J. Indeed, the integrals that arise when taking the weak-coupling limit, such as 

t t 

j drF z s T/x2 = lj^drj dug^P^a^e-^^ 2 , (20) 

from Eq. f[T9l become extremely singular as A — > 0. To get a feel for this singularity, it 
is useful to consider the commutator 

lim [s t/X 2, 4 /A2 ] = 2ir\g(uj A )\ 2 5(t - s), (21) 

which suggests that the field operators are delta-correlated in time. Furthermore, the 
quadrature operators (s t — s[)/i commute at different times, which suggests that one 
should interpret the field operators as white noise, viewed heuristically as arising from 
the vacuum fluctuations in the probe field. This analysis can be made mathematically 
rigorous, and the interested reader is encouraged to consult, for example, Refs. [7j [27]. 

Unfortunately, the presence of white-noise drive terms in Eq. ( fT9l) means that the 
integrals do not converge in the same way that they would were the integrand a smooth 
function. One must therefore employ the methods of stochastic calculus [32J to evaluate 
the 5^0 limit. By analogy to the procedure adopted in classical stochastic, the 
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extreme singularity of the the quantum white noises s t and sj, is alleviated by defining 
the quantum Brownian motion operator 

S t = f Kdr (22) 



which enables one to define the the quantum Ito integral 

f X T dS T = lim £ X tk (S tk+1 - S tk ) (23) 

where \tk+i — tk\ =8 and X is an arbitrary atomic operator. The dSj and dSt are 
delta-correlated noise operators derived from the quantum Brownian motion and satisfy 
the quantum Ito rules: dS t dS\ = dt and dS\dS t = dS 2 = (dS}) 2 = 0. It is a byproduct 
of the white-noise singularity that the quantum Ito integral does not obey the regular 
chain rule of differential calculus, but rather the Ito rule 

_ df(X t ,t) df(X t ,t) 1 &f(x t ,t) 2 
df - —Qx^dX t + —g^dt + - Qx? dX t , (24) 

where f(X t ,t) is any well-defined function of X t and t. This modification of the 
chain rule is the price one must pay to obtain noise increments that are statistically 
independent from the state of the system E j* X T dS T = 0. Sometimes however, it is 
desirable to retain the regular chain rule of calculus, 

<f = S J§f> dX , + S J^ dt , (25, 

such as will be the case when applying the methods of differential geometry to stochastic 
systems [as will be considered in Sec. (13. 2ft ], The regular chain rule of differential calculus 
can be recovered by defining the stochastic integral as 

I Xt o dS T = lim ]T \ (X tk+1 - X tk ) {S tk+1 -S tk ), (26) 

which is called a Stratonovich integral. However, E J * X T o dS T ^ 0, making it generally 
less desirable for analysis than its Ito counterpart, Eq. (|2"3"|) . 

With a well-defined quantum stochastic calculus given by Eqs. fl23|) and (126]) . it is 
possible to obtain the first- and second-pass atom-field interactions. The result of the 
derivation is the following quantum stochastic differential equations (QSDE) [22] 

dU t {1) = ^/WxF^dSt - dS t ) - X -m x Fldt - ^H t dtj (27) 

dlli 2) = iyiy/W 2 F y {dS t + dS\) - X -m 2 F 2 dt - l jHtdt^ (28) 

where mi = 27r|Ajg(u;A)| 2 . Here Ht is an arbitrary Hamiltonian applied to the atoms 
that is not involved with their interaction with the probe field. 
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2.3.2. Combined Weak-Coupling Propagator. To obtain a single weak-coupling limit 
for the double-pass interaction, the separate equations of motion for the two propagators 
must be combined into a single weak-couping limit. To obtain such a coarse-grained 
propagator, the two single-pass evolutions dU^ — d t U^ and dU^ = b t U^ can be 
concatenated via the expanded differential 

dU t+ 2St = Ut + ( a + b + bo) U t (29) 



such that the combined propagator satisfies d\J t = Ut+2&t — Ut in the limit where 25t is 
taken to be infinitesimal|§|. After evaluating the combined evolution for the propagators 
in Eqs. (j2"7|) and (1281) in light of the quantum Ito rules, the single weak-coupling limit 
propagator is found to be [22] 

1 , »■) , i 



dU t 



VKMF y F z dt - -MF^dt - -KFHt - -Hdt (30) 



y z 2 2 y h 

t _ JC A i „- A /zFl? Met 



+VMF z (dSj - dS t ) + WKF y (dS] + dS t 



where M = m>i/2 and K = m 2 /2. Finally, this result can be adapted to the standard 
notation of quantum stochastic calculus by rewriting the double-pass propagator in the 
general form 



dU, 



LdS\ - VdS t - -Dtdt - iH c dt 



U t , (31) 



with the assignments L = VMF Z + i\fKF y and H c = H — \ / KM(F z F y + F y F z )/2 

[221 EH]- 

The stochastic propagator allows one to compute the time-evolution of any atomic 
or field observable. For example, the Heisenberg-picture operator for any atomic 
observable X (time-independent in the Schrodinger picture) is given by the quantum 
flow jt(X) = U\{X ® I)Ut whose time evolution satisfies 

^(X)=j 4 ([X,L])^+j 4 ([Lt,X])^+j t (£(X))rft (32) 

where C denotes the familiar Lindblad generator for the operator X, 

C(X) = i[H, X] + DXL - -{i^LX + XL r L). (33) 
3. Continuous Observation and Quantum Filtering 

In addition to the stochastic propagator U t , it is also useful to consider the effect of 
performing measurements on the twice forward-scattered laser field, as illustrated in 
Fig. ([I]). The task is then to construct a best-estimate of atomic observables conditioned 
on the measurement data, which can be accomplished using the techniques of quantum 
filtering theory [TJ [27]. Of course, the form of any such filter depends on how the 
forward scattered probe field is detected. From the form of the first-pass interaction 
Hamiltonian \ it is evident that the ^-component of the atomic spin couples to 

§ Care must be taken to account for rescaling time from 2St — > dt, as described in Ref. [55] 



Coherent Feedback in Collective Spin Systems 



11 



dynamics generated by the field operator s Zjt = i(s\ — s t ). The effect of such a 
coupling would then be observed by measuring the orthogonal quadrature, indicating 
that information can be gained about the ^-component of the atomic spin by measuring 
the polarization observable z t = s y t [331 122] . The observation process considered here 
is therefore described by 

Z t = u}(Sl + S t )U t . (34) 

whose time evolution satisfies the quantum Ito equation 

dZ t = ] t (L + D)dt + dS\ + dS t . (35) 

In the language of classical control theory, Eqs. ( 1321) and ( |35l) provide the system- 
observation pair. The observations process carries information about the state of the 
system, j t (L + U), albeit corrupted by the quantum white noise term dS} + dS t . 

3.1. The Quantum Filtering Equation 

Quantum filtering theory is tasked with picking out the relevant information from the 
observations and combining it with the model of the dynamics in order to estimate 
system observables. Mathematically, the filter is given by the conditional expectation 

7r t [X] = E[ Jt (X)\Z m }, (36) 

where Z(o,t) is the measurement record up to time t. It can be shown that the conditional 
expectation in Eq. (136]) is adapted to the measurement process and that it is possible 
to express the evolution of the quantum filter as an Ito equation 

dir t [X] = a t (X)dt + b t {X)dZ t , (37) 

where at(X) and bt{X) are functions determined by enforcing 7r t to satisfy the properties 
of a conditional expectation, namely 

E[YE[X\Z m }} = E[YX] (38) 

for any Y that is a function of the measurement process. Doing so yields the quantum 
filtering equation for the double-pass atomic system [22] 

dn t [X] = n t [C[X]]dt + (ir t pX + XL] - n t [P + L]vr t [X]) (39) 

x (dZ t - ir t [L + ti)dtj . 

The quantum filter provides the best least-squares estimate of the expectation value of 
any atomic operator X given the measurement data. This suggests that it should be 
possible to find a density operator p t such that vr t [X] = tr[Xp t ] for any observable X. 
Such an operator can indeed be found [3H [7], and it satisfies the adjoint form of the 
quantum filter 



dp t = - l j[Hu Pt]dt + iVKM[F y , {F z , p t }]dt (40) 



+ MV[F z }p t dt + KV[F y ]p t dt + VMM [F z ]p t + iVK[F y , p t ) dW t 
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where the innovations process 

dW t = dZ t - 2\fMtr[F z p t ]dt (41) 
is a Wiener process, i.e. E[dWt] = 0, dW? = dt, and the superoperators are defined as 

V[F k ]p t = F k p t Fl - \PlF k p t - l -p t FlF k (42) 

M [F z ]p t = F z p t + pA - 2ti[F z p t ]p t (43) 
{F z ,p t } =F z p t + p t F z . (44) 

In situations where the conditional density operator remains pure under the evolution 
of the adjoint filter, as is the case for evolution under Eq. (j4*0l provided that the initial 
state is pure, p t can be decomposed as p t = iV'tKV'tl and shown to satisfy the stochastic 
Schrodinger equation (SSE) [21] 

i\fKMF y [F z + (F z ) t ) - f(F z - {F z ) t f - ^3 - ^ 

+ (VM(F Z - (F z ) t ) + iVKF y ^) \^ t )dW t . 
where the superoperator M is defined as the generator the dynamics. 

3.1.1. Amplification of y-Axis Rotations. The first term in Eq. ( )45l) suggests the 
appearance of a Hamiltonian-like term that drives rotations about the y-axis in a manner 
that depends upon the current z-component of the atomic spin. Now consider a system 
Hamiltonian that involves a y-axis rotation 

H t = -hu t F y (46) 

where u t is a time-dependent rotation rate about the atomic y-axis. The filtering 
equation then becomes 




M ~ ,~ KF 



2 



d\^ t ) = \ iu t F y + iVKMF y (F z + (F z ) t ) - — (F z - (F z ) t ) 2 - ) \^ t )dt 



+ (VM(F Z - (F z ) t ) + tVKF^j \ip t )dW t . (47) 



Since the sign is the same for both the drive and feedback terms (first and second terms 
in the filtering equation), the resulting feedback is positive and amplifies the effect of 
the Hamiltonian Eq. (|4"6"|) . The generator of the feedback dynamics is not, however, 
Hermitian because the process is not reversible. The feedback on F y is dependent upon 
F z via the interaction of the system with a reservoir. As in the theory of cascaded 
quantum systems [35], both the first and second-pass interactions are one-way due to 
the immensity of the bath. 

Figure (J2j) illustrates the amplification process just described by comparing the 
dynamics for a spin system with F = 100, uo t = 2.0 and M = K = 1.7 (with frequencies 
in units of 1/tfmai) to one in which the second-pass interaction strength is set to zero 
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Figure 2. (color online) Comparison between the atomic dynamics with and without 
positive feedback amplification enabled. The measurement trajectories depict the filter 
TT t (F z ) = E[F z |Z C0>t) ] = {ip t \K\il>t) evolved under Eq. g7|) with F = 100, w t = 2.0 and 
M = K = 1.7 (in units of 1/tfinai)- The double-pass filter shows an amplification 
of the single-pass dynamics, which were obtained by setting K — for the sake of 
comparison. 



K = 0. The filtered value iTt(F z ) is plotted for both cases, and the nonlinear effects of 
the positive feedback are evident from the trajectory. In order to obtain a meaningful 
comparison between the evolutions with and without feedback amplification enabled, 
the filter Eq. (H71) was propagated using the same noise realization for the innovations 
process dW t , which is possible in a computer simulation by using the same random seed 
for both cases. Of course, the amplification process will also affect the "quantum noise" 
in the atomic system, due to the Heisenberg uncertainty of the initial spin coherent state. 
The benefit of amplifying a signal Eq. ( ]46i) relative to also amplifying the coherent state 
projection noise is analyzed in the remaining sections. 

3.2. Reduced Dimensional Filters 

While the quantum filter Eq. ( 1451) provides a complete description of the conditional 
dynamics of the double-pass atomic system subject both to the coherent feedback 
mediated by field and to continuous observation of the field, it is somewhat impractical to 
utilize in practice. Typical experiments involve atomic spin systems with N ^> 1 atoms. 
Despite that the effective atomic Hilbert space, the completely symmetric representation 
described in Sec. (J2]), grows only linearly with N, the corresponding operators and kets 
in Eq. ([4"5]) are still too large to be analyzed numerically on a typical high-end computer 
unless F is on the order of a few hundred or less. 

Preferably, it would be possible to parameterize the atomic Hibert space even 
further, ideally in a way where the number of parameters does not increase with the 
number of atoms in the atomic sample. Such a parameterization is provided by the set 
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of Gaussian states [351 E3 EH] that arise when it is possible to truncate the Holstein- 
Primakoff expansion [39] of the atomic spin state \ip t ) to first order (as described below). 
Such states arise naturally when an atomic system evolves from an initial spin coherent 
state, such as the spin-polarized states obtained by optical pumping [40J. Then, an 
effective filtering equation can be derived for the small number of parameters that 
characterize the state. 



3.2.1. The Family of Gaussian States and Restricted Dynamics. Such a reduced- 
dimensional filter can be obtained using the two-parameter family of Gaussian states 

\9 t ,^)=R y (9 t )S^ t )\F,+F) x , (48) 

where \F,+F) X is the spin coherent state pointing along +x, S(£t) is a spin squeezing 
operator [UJ 

%) = exp (-2i£ t {F z F y + F y F z )) (49) 
with squeezing parameter £ f and R y {0t) is a rotation about the y-axis 

R y {9 t ) = exp (-i0 t F y ) (50) 

by angle 9 t . This parameterization is, of course, limited; it can only capture dynamics 
that involve rotations about the atomic y-axis and the squeezing that is generated along 
the z-axis that results from measuring F z . The rotation via R y (9 t ) then accounts for 
both the random evolution due to the measurement and the feedback term that drives 

3.2.2. First-Order Holstein-Primakov States. For large F, spin-polarized states are 
extremely well described by the Holstein-Primakoff approximation to lowest order [39J 

F +>x « V2Fa (51) 

F_ x ss v / 2Fa t (52) 

F x P3 F, (53) 

where F± )X = F y ± iF z , and d^,a are bosonic creation and annihilation operators. 
The state can then be written as \F,+F) X = |0), which is the vacuum in the 
Holstein-Primakoff representation. When working with this parameterization, one often 
encounters inner products of the form 

X (F, +F\S&y g (F x , F y , F a )4(0 t ) + /(i^ F y , FjRyiOJSfo) \F, +F) X , (54) 

where g and / are polynomial functions of their arguments. Since R y {9 t ) is a rotation, 
it is possible to evaluate 

R y (9 t yf(F x , F yi F z )R y (9 t ) = f(F>, F y , F z '), (55) 

exactly using 

F x = R y (9 t yP x R y (9 t ) = F x cos 9 t + F z sin 9 t 

P; = R y (9 t )lF y R y (9 t ) = F y (56) 
F' z = R y {9 t yP z R y {9 t ) = F z cos 9 t - F x sin 9 t . 
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Thus it only remains to evaluate matrix elements of the form 

X (F, +F\S(^g(F x , F y , F z )f(F^, F y1 F')S^ t )\F, +F) X (57) 

which is possible for small £ 4 , i.e., for "squeezed vacuum" in the preferred basis. Under 
this approximation 

S(Z t yP x S(Z t ) = F (58) 



S(0F y S(Z t ) = ^_ e +4 ^ (a + fit) (59) 
SfayFJSfa) = ^-e" 4 ^(a - fit). (60) 

3.2.3. Quantum Projection Filter. With the Gaussian family of states defined, an 
approximate quantum filtering equation confined to the family can be derived. Our 
approach utilizes the technique of projection filtering [4"2l IT?] . Abstractly, the quantum 
projection filter operates by orthogonally projecting Eq. (|4"5T) onto the sub-manifold of 
states defined by Eq. f|48l at each point in time. In doing so, the objective is to obtain 
an equation of motion 

d\^ t ,e t ) = v^ t + ve t d6 t . (61) 

in terms of the parameters of the family, namely 9 t and £ t , rather than an equation of 
motion for the full ket \ipt). 

Given a family of states characterized by the set of parameters x±, x 2 , ■ ■ ■ , x n and 
with states in the family \xi,x 2 , ■ ■ ■ ,x n ), the tangent space for that family is spanned 
by the vectors [12] 

d\xi,x 2 , ■ ■ ■ ,x n ) 

Vi = ^ ' (62) 

which describe how states within the family \x%, x 2 , . . . , x n ) are modified under 
differential changes in the values of the parameters. Of course the action of the generator 
jV[|xi, x 2 , ■ ■ ■ , x n )] such as that defined by Eq. fH5|) does not necessarily preserve states 
in the family. But, projecting Af onto the tangent space 

rp (Vi,M'[\Xl,X2,...,X n )]) , 

where the inner product is given by the standard Hilbert space inner product, yields 
the filter that is the closest approximation to Af that preserves states within the 
parameterized family. For the family of Gaussian states defined by, the tangent vectors 
are found to be 



v 9t = d -^- = -iF y R y (9 t )S(Ct)\F, +F) X (64) 
«6 = = R y (e t )Sm-^(F z F y + F y F z ))\F,+F) x . (65) 
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with the inner products required for normalization of unit vectors given by 

Fe SF£ t 



(V0t,V0t) 



8F 2 
0. 



(66) 



The last line of Eq. (1661) indicates that the tangent vectors are orthogonal as required. 

3.2-4- Orthogonal Projection of Double-pass Filter. To construct the projection filter 
onto Gaussian states for the double-pass quantum system in Fig. (pQ), all that remains 
is to implement the projections. Unfortunately, this cannot be accomplished using Eq. 

directly as the quantum filter in this form is an Ito equation, whose chain rule is 
incompatible with the conventional methods of differential geometry [7]. Fortunately, 
Stratonovich stochastic integrals obey the standard chain rule and are thus amenable 
to projection filtering methods. Transforming Eq. (1451) to Stratonovich form can be 
performed using the approach described in, for example, Ref. [22] and gives 

\fKM - 



-iuj t F Y - M 



(F z - (F z ) t y - (AF Z ) 



+2iVKM{F z ) t F y + iVKM{F z F y ) t 
+ (VM(F z -(F z ) t )+tVKF y 



\ij}) t dt 
odW t , 



(67) 



where (AF z ) t = (F z ) — (F z ) 2 . The projection filter is then obtained by comparing Eqs. 
fl61~]) and fl63]) using Eq. fIBTl) and the general forms for d£ t and d6 t 



d6 t 



1 



8F 2 



>*,d|Vt>Mt]> 
(vt t ,d\&)[&A]), 



(68) 
(69) 



where d\ipt)[£,t, &t] is the evolution of |£ t , 8 t ) under the Stratonovich filter of Eq. (ISTj) . 
After converting the resulting expressions back to an Ito equation from its Stratonovich 
form, the resulting two-parameter filter on Gaussian states is given by 

M 



dd t 



0J t 



dt 



(70) 



Me" 8ff ' cos6 t + yfE 



dW t 



M 



e" 8i ^ cos 2 9 t dt, 



where the innovations are also expressed within the Gaussian state approximation of 
the filtered expectation value (F z )t = —Fsm9 t to give 

dW t = dZ t - 2VM(F z ) t dt 

= dZ t + 2FVM sm 6 t dt. (71) 
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Measurement Time (units of \IM) Measurement Time (units of MM) 

Figure 3. (color online) Comparison between the projection and full quantum filtering 
equations for the measurement trajectories in Fig. ([2]). In plot (a), 7r ( (F z ) = E[F z |Z t ] 
is computed for a simulated measurement trajectory using both the projection filter 
Eq. (|70|) and the full filter Eq. ((47)) under the same conditions as used in Fig. (T5J). 
The two filters given nearly indistinguishable resulsts. Plot (b) shows the relative 
error |(7Tj 6xact ' ) (F z ) — 7T( Pro ^)/7rj Cxact ^(F z )| for both the double-pass and single-pass 
evolutions. 

Figure ([3]) illustrates the accuracy of the projection filter by comparing the 
conditional expectation value 7r t J '^(F Z ) = E t [F z |Z(o,t)] computing from Eq. ( 1701) with 
the conditional expectation 7r ^ cxact ' ) obtained from Eq. (|4"TI) under the same conditions 
used to generate Fig. (TSJ). The relative error in the projection filter is found to be below 
1%, except for regions where the filter crosses zero (which accentuates the relative error). 
One would expect the projection filter to perform increasingly better with larger values 
of F, but the comparison in the figure is limited to values for which Eq. ( H7|) can be 
evaluated in practice. 

4. Quantum Parameter Estimation and Atomic Magnetometry 

Since the effective dynamics of the double-pass atomic system in Fig. (TjQ) are seen to 
amplify rotations about the y-axis, it raises the question whether such a configuration 
provides an improved capability to detect external influences that drive y-axis rotations 
in the atoms. In atomic magnetometry [121 SUES], for example, a y-axis magnetic field 
B = B y is typically detected by preparing an atomic sample into the x-polarized initial 
state \F,+F) X , allowing it to undergo Larmor precession until some time t and then 
inferring the value of B by measuring the ^-component of the atomic spin to determine 
the amount of rotation. For very precise measurements, uncertainty 5B in the estimated 
value B of the field is dominated by quantum fluctuations in the observations performed 
on the atomic sample [15J H6j H7] . 

Magnetometry fits naturally into the coherent feedback plus continuous 
measurement structure developed in Sections ([2]) and (J3j). The atomic sample couples 
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to the magnetic field via the Zeeman Hamiltonian 

H = -hjBF y , (72) 

where 7 is the atomic gyromagnetic ratio, which has precisely the structure considered 
in the derivation of the Gaussian state projection filter with uj t = 7-B in Eq. (146|) . 
Measurements of the ^-component of the spin reveals information about Larmor 
precession, which should then in-turn enable estimation of B by devising a suitable 
estimation procedure [as addressed in Section (14.31) ]. 



4-1. Lower-Bounds in Quantum Parameter Estimation 

The first problem at hand, however, is to ascertain whether the positive feedback 
amplification of Eq. (!72|) offers a potential advantage for determining the value of B 
as compared to an estimation procedure that relies on only a single-pass measurement 
configuration [38]. At first brush, it would seem that amplification of the Larmor 
dynamics such as observed in Fig. (j2J) should make the influence of the field easier 
to observe, especially for small values \B\ <C 1. Amplification of Eq. (!72|) is, however, 
not the entire story because an amplification of F z occurs even in the absence of an 
applied field. A magnetic field estimator must be able to distinguish the amplification 
of the magnetic signal from amplification of the coherent state spin projection noise 
[3H S3 HE]. 

The potential for positive feedback amplification to improve the parameter 
estimation uncertainty can be analyzed via the quantum Cramer-Rao inequality 
[4"5 j |4~9 | \%6 \ [50] . which places an information-theoretic lower bound on the units-corrected 
mean-square deviation of any estimator B t from the true value of the field B, 

1/2 

SB t = ( I b) ) (73) 

\\\d(B)/dB\ 1 1 

in terms of the dynamics 



5B t > tr 



4MB) 



dpt(B) \ 
dB ) 



(74) 



The behavior of the estimator uncertainty depends on the characteristics of the quantum 
states used to compute the expectation value in Eq. (IT4"]) and also the nature of the 
induced dynamics [51]. If one considers only Hamiltonian evolution under Eq. (!72j) and 
does not permit quantum entanglement between the different atoms in the probe, the 
Cramer- Rao inequality then yields the so-called shotnoise uncertainty [US [38] 

5^ = -^=, (75) 

whose characteristic l/V~N scaling is a byproduct of the projection noise (AF Z ) = 

for a spin coherent state [HI SB]- If, however, entangled states are allowed to 
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evolve under Eq. (I72"j) . then the Cramer- Rao inequality yields the so-called Heisenberg 
uncertainty [HI EH EH [52] . 

SM HL) = — *— , (76) 

which offers a quadratic improvement over shotnoise as a function of the number of 
atoms. Furthermore, if one goes beyond the structure of Eq. ( 1721) by implementing 
dynamics that are non-linear in the collective spin operator, it is possible to surpass the 
1/N scaling [51] as a consequence of multi-body effects [5Tj [53]. 

4-2. Numerical Evaluation of Cramer-Rao Bound 

A lower-bound on the magnetic field estimation uncertainty possible under the dynamics 
generated by the double-pass configuration in Fig. ([1]) can be obtained by integrating 
the quantum filter in Eq. (pf5l) and computing a finite difference approximation to the 
operator derivative dpt(B)/dB in Eq. (17^|) . Such a finite difference approximation can 
be constructed by co-evolving three trajectories, f>t(B), p t (B + 5), and p t (B — 5), and 
calculating 

' Mb + s) - pt{B - sy 




Pt(B) 



25 



(77) 



for 5 <C 1. The estimator uncertainty bound 5B t obtained in this manner is 
conditioned on the random measurement realization driving the innovations process 
dWt for any single integration of the quantum filter. Averaging Eq. ( 1771) over many such 
measurement realizations then provides the unconditional bound corresponding to the 
double-pass magnetic field estimator. 

A thorough numerical investigation of the behavior of 5B t computed from the 
Cramer-Rao inequality Eq. ( 17^1) evaluated using Eq. ( 1771) reveals that the estimator 
performance depends very highly on the values of the coupling parameters M and K, 
and that the best values of M and K depend on the spin size F [2TJ [22]. This is 
not surprising, as even in a classical amplifier, achieving increased signal to noise often 
requires selecting an optimal value for the amplifier gain. The current situation is no 
different. By conducting large sets of simulations, we identified empirically that scaling 
the values of M and K according to the functional form [21] 



M — K — — %- , (78) 

where c and a are constants, leads to a power-law scaling for the uncertainty bound 
5B tfiual ~ 1/N k . As demonstrated by the data points in Fig. (HH), it is possible to 
achieve 1/N scaling (to within a small prefactor offset) by setting a ~ 0.77. A linear 
fit to the Cramer- Rao data (on a log- log scale) gives SB t&ii&l ~ ^-o.98^ w j 1 j c i 1 reproduces 
the Heisenberg scaling to within the precision of the fit. The data in Fig. (jib) shows 
the slope of a linear fit of \og 10 5B tfin!ll to log 10 F (i.e., a slope of k = — 1 corresponds 
to the Heisenberg uncertainty scaling) for different values of a in Eq. ([TBI) . This plot 
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Figure 4. (color online) Plot of a lower bound on the estimator uncertainty for atomic 
magnetometry performed according to the measurement configuration in Fig. (TT]) as 
determined from a numerical evaluation of the quantum Cramer-Rao inequality. Plot 
(a) demonstrates that it should be possible to achieve 1 /F scaling of the magnetometer 
sensitivity by adjusting the measurement parameters M and K as a function F 
according to Eq. (f78|) with a — 0.77. A fit to the data points in plot (a) gives a 
scaling law of i^-°- 98 ; which is consistent with 1/F to within the accuracy of the fit. 
Plot (b) illustrates that is possible to surpass the Heisenberg uncertainty scaling using 
different for different values of a. 



suggests that it should be possible to surpass the Heisenberg "limit" as well via a double- 
pass measurement, however, these calculations provide no guarantee that the Gaussian 
state approximation continues to hold for larger values of M and K. In any given 
measurement configuration, it is essential to ascertain whether or not the Gaussian 
family provides a sufficient approximation to the actual atomic state. 

Extrapolating the values of M and K that would be required to achieve a sensitivity 
<5i?t flnal corresponding to the Heisenberg uncertainty scaling Eq. ( 1761) with N ~ 10 8 
(which is a conservative value for typical experiments) implies that one would only 
require Mt final ~ 10~ 7 , which suggests that 1/N scaling can be achieved by an extremely 
weak measurement, with no appreciable generation of conditional spin-squeezing [21J. 
This suggests that it may be much easier than previously believed to improve the 
sensitivity of an atomic magnetometer in practice — scalings beyond 1/yN are made 
possible through the use of amplified dynamics instead of the preparation of an entangled 
state of the atoms. In the case of small M and K, the Gaussian family should provide 
a good representation of the atomic state. 

4-3. Quantum Particle Filter 

The Cramer-Rao inequality provides an ideal lower-bound on the sensitivity of 
a parameter estimator, but does not generally yield a constructive procedure for 
implementing parameter estimation. It is, however, possible to obtain a magnetic 
field estimation procedure for the double-pass measurement using a technique referred 
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to as quantum particle filtering [22] • The concept of quantum particle filtering was 
developed in Ref. [54]. It operates by embedding a classical probability distribution 
for the parameter to be estimated into a quantum probability space [HU [22]. This is 
accomplished by promoting the classical parameter to be estimated — in this case B — 
to a diagonal operator 

B ^ B = J B\B)(B\dB eH B , (79) 

where H,b is an auxiliary Hilbert space for the parameter with basis states satisfying 
B\B) = B\B) and (B\B') = 5{B — B'). All atomic operators and states, which are 
associated with the atomic Hilbert space Ha, act as the identity on Hb, e.g. F z ^ I<S>F Z . 
The only operator which joins the two spaces is the parameter coupling Hamiltonian, 

H ^ -HjB ® F y . (80) 

The standard methods of quantum filtering theory are then employed to obtain the 
best least-squares estimate of the parameter via the conditional expectation 

B t = n t (B) = E[ Jt (B)\Z m } (81) 

following the approach detailed in Section (J3]). Since B corresponds to a classical 
parameter, the marginal density matrix tr^ A [p t ] should be diagonal in the basis of B so 
that it represents a classical probability distribution. Therefore the conditional density 
operator propagated by the adjoint form of the quantum filtering equation 

Pt = J Pt{B)\B){B\®p[ B) dB (82) 

where Pt{B) = P(B\Zm t )) is precisely the conditional probability density for B. The 
density function describes a continuous variable, but in practice is discretized using a 
weighted set of point masses or particles: 

N 

Pt (B)^j2pt )s ( B - B ^ ( 83 ) 

i=l 

an approximation which can be made arbitrarily accurate in the limit of iV — > oo. Under 
this discretization, the conditional density operator in Eq. ( |82i) takes the form 

JV 

pf^Y.^WiB^p^, (84) 

i=l 

where each of the iV triples {pf \ Bi, Pt } is called a quantum particle, giving rise to 
the name of the parameter estimation procedure. The parameter estimate can then be 
obtained from the approximate density in Eq. ( |83i) either by taking the most probable 
B value, i.e., the Bi that maximizes pf , but generally by calculating the expected value 
of B 

B t = n t {B) (85) 

N 



i=l 
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and the parameter estimation uncertainty is obtained from the variance of the 
conditional density 



N 



6B* = (B?) - B\ 



(86) 



i=i 



4-3.1. Gaussian State Particle Filter. The quantum particle filter for the double- 
pass system within the Gaussian family restriction is found by substituting p t into 
the extended double-pass filter and then projecting onto the manifold Eq. (HHj) . The 
resulting particle filter is given by 

N 

i=i 



dp 



(0 



2FVM 



sin 9\ 



(BO 



<0) 



sin 6 



pfdW t 



87) 



dd 



7^ - ^ ex P (- 16i? & (Bi) ) sin20 t (Bi) + 2FVKMsm6 { t Bt) 



dt 



M exp (-8F& (Bi) ) cos0 ; 



j(Bi) 



^ = ^exp(-8FeP ) )cos^f^t, 



with the modified innovations process 



A? 



dWi = + 2FVM^pf } sin^ (Bl) dt. 



(88) 



i=i 



Evidently, each particle's quantum state p t evolves under the projection filter 
corresponding to the known field B = Bi, while the innovations process dW t serves 
to couple the different particles. The form of this filter, which involves several terms 
that are nonlinear in the variables 9 t and £ 4 suggest that, even if the state is Gaussian 
and the estimator density Pt{B) is of a Gaussian form, it is still unlikely that Kalman 
filtering methods [38] can be applied as a result of the nonlinear dynamics. 



4-3.2. Parameter Estimation Simulations. We performed simulations of the Gaussian 
projection particle filter magnetic field estimator Eq. (187]) to assess its performance 
relative to the Cramer- Rao predictions in Fig. (j4]). As is often the case when analyzing 
detection limits, we took the true value of the field to be B = 0, and evolved a statistical 
ensemble of measurement realizations under Eq. <\87L For our simulations, we utilized 



S9) 



N p = 10 4 particles and took the prior distribution over B to be uniform 



(i) 

Po 



with the particle values Bi placed on a uniform grid with domain ±y D centered around 
B = 0. That is, the particles have the values 



B 



j-v^D, -VD + dB,...,VD-dB, Vd} 



(90) 
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Figure 5. (color online) Numerical simulation of the Gaussian projection particle filter 
for estimating the value of an external magnetic field that uses coherent amplification 
of the magnetic dynamics to achieve a sensitivity that scales as 1/F without the need 
for spin squeezing. 



with dB = 2\^D/N p . For our dynamics, we considered 7 = 1 with the initial quantum 
state given by \9 ,£, ) with 9 = and £ = 0, i.e., the x-polarized spin coherent state. 
The coupling parameters, M and K, were scaled according to Eq. ( 1751) using c = 0.5888 
and a = 0.77, whose values were chosen such that the estimator uncertainty should 
scale as 1/F with the size of the spin system. 

The results of the simulations are given in Figure (jSJ), which plots the estimator 
uncertainty SB obtained from Eq. ( 1861) as a function of F. The data points indicate the 
average value of SB over the statistical ensemble of 10,000 measurement realizations, and 
the error bars indicate the standard deviation of SB over that ensemble. For reference, 
the shotnoise and Heisenberg uncertainty scalings are also plotted on the axes, making 
it apparent that the projection particle filter clearly outperforms shotnoise. A numerical 
fit of the projection particle filter scaling gives SB ~ ^-o.98^ ^jgjj j^g fj^g same power 
law slope as the Heisenberg scaling to within the precision of the fit. The absolute 
magnitude of the parameter estimation uncertainty SB exhibits, however, a displacement 
from that of the Heisenberg uncertainty limit. This is not surprising, since the Crame- 
Rao inequality provides a theoretical lower bound, and there is no guarantee that an 
estimator (even an optimal estimator) is capable of achieving that bound. As expected, 
based on the analysis of Section fj4.2j) . the atomic state not only remains Gaussian, but 
also undergoes extremely little squeezing. For the measurement realizations considered 
here and scaling under Eq. (1781 . typical values of £ were found be £ ~ 1.2 x 10~ 3 for 
F = 100 and £ ~ 2.0 x 10~ 4 for F = 1000. Since £ = corresponds to a spin coherent 
state, this constitutes an improvement over shotnoise without appreciable generation 
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Figure 6. (color online) The conditional probability density produced by the 
projection particle filter exhibits a bias unless a sufficient number of particles is 
employed and the prior distribution is sufficiently wide. The figure demonstrates 
reduction of bias and convergence of the estimator as the width of the prior distribution 
is increased. In each case, the prior distribution is taken to be uniform over the range 
of S-values given by ±\^D. The filter is seen to converge (for this measurement 
realization) as D is increased from 1 x I0 3 to I x I0 5 . The solid line is a fit of the 
evolved conditional probability density to the converged normal distribution. 

of entanglement between atoms. Extrapolating our results to experimentally viable 
values of F ~ 10 8 , or more, suggests that coherent feedback amplification of the Larmor 
dynamics enables 1/F uncertainty scaling in parameter estimation without the need for 
spin squeezing, making the approach attractive for practical application. 

Several important factors remain be to be solved, however, before the procedure 
developed here could be successfully incorporated into an actual laboratory setting. As 
in any situation where a continuous probability distribution is discretized, as required for 
numerical analysis, there is the potential to introduce bias into the estimator outcome, 
even if the underlying analytic methods are bias-free [541 [55] . During the stochastic 
evolution of the particle filter, especially the early transient period of the evolution, it 
is possible for the distribution mean to vary significantly. Noise realizations that couple 
to particles which lie far away from the true value of B can develop a bias if the tails of 
the distribution exceed the domain to any significant degree. Such a bias might offset B 
from the true value of B by an amount that exceeds the estimator uncertainty 5B, and 
thus prevent the estimator from achieving its expected uncertainty. In an experiment, 
such a bias would be diagnosed by observing that the sample variance of the estimator 
E[(B - B) 2 ] exceeds that predicted by Eq. (|86l) . 

We have observed such a bias in our simulations, and have found that care must be 
taken to minimize the introduction of bias by using a sufficiently wide prior distribution 
on B. Figure ([6]) demonstrates the reduction of estimator bias as the width of the prior 
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distribution is increased from D = 1 x 10 3 to D = 1 x 10 5 . With increasing D, the mean 
of the estimator distribution converges toward the Gaussian distribution illustrated 
by the solid line in Fig. (jHJ). We have also found that the value of D required for 
such convergence exhibits some variability with respect to the particular measurement 
realization. This is not unexpected, as the degree to which the filter must incorporate 
particles that lie near the domain boundaries depends upon the stochastic propagation, 
and different noise realizations cause the filter to explore 5-space differently. 

Increasing the domain by taking D ^> 1 does not come without a cost. For a fixed 
number of particles, the filter grid size dB = 2\^D/N p increases with D. For a fixed 
value of F, a smaller fraction of the particles have non-negligable weight at the end of 
the filter evolution as D is increased. This reduced resolution is also evident in evident 
in Fig. ([6])- the larger domain sizes coincide with reduced resolution of p(B\Z^). The 
number of particles must be chosen such that there is adequate resolution p(B\Z(p tt )) to 
determine its mean and variance with confidence. As the number of atoms increases, 
F the width of p(B\Zf 0)t \) decreases, suggesting that the number of particles must also 
be increased to maintain adequate resolution. As N p increases, it becomes increasingly 
time-consuming to integrate the filtering equations, and this could prove limiting in an 
actual magnetometry application. Fortunately, there are a number of potential methods 
one might use to improve the efficiency of the filter by reducing the effective number 
of particles. For example, particles whose weight drops below a given threshold could 
be removed from the integration, or a more sophisticated basis than the delta function 
representation Eq. (183]) could be used to resolve p(B\Z^). While a full exploration of 
the numerical methods suitable for optimizing the performance is beyond the scope of 
this work, such practical considerations are likely to play an essential role in deploying 
particle filters in a laboratory experiment. 

5. Conclusions 

We have presented a method for quantum parameter estimation, especially suited 
toward atomic magnetometry, based on the amplification of Hamiltonian dynamics via 
coherent positive feedback. We also explored the behavior of our feedback system using 
information-theoretic parameter estimation bounds, via the Cramer-Rao inequality. To 
render our approach more suitable for incorporation into an actual laboratory setting, 
we employed methods from quantum filtering theory and differential geometry to obtain 
a simplified estimator based on a two-parameter family of Gaussian atomic spin states. 
Finally, we demonstrated that it should be possible for an atomic magnetometer to 
achieve a sensitivity scaling comparable to the Heisenberg scaling without the need 
for entanglement. Our results suggest that effective nonlinear dynamics, such as those 
obtained here through positive coherent feedback, are likely to play as important a role 
in the field of quantum parameter estimation as classical amplifiers play in classical 
metrology. 
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